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ABSTRACT 

We present timing observations of four millisecond pulsars, using data obtained 
over three years at the ATNF Parkes and NRAL Jodrell Bank radio telescopes. As- 
trometric, spin, and binary parameters are updated, and substantially improved for 
three pulsars, PSRs J0613-0200, J1045-4509 and J1643-1224. We have measured 
the time variation of the projected semi-major axis of the PSR J0437— 4715 orbit 
due to its proper motion, and use it to constrain the inclination of the orbit and the 
mass of the companion. Some evidence is found for changes in the dispersion measures 
of PSRs J1045— 4509 and J1643— 1224. Limits are placed on the existence of plane- 
tary mass companions, ruling out companions with masses and orbits similar to the 
terrestrial planets of the solar system for eight pulsars. 
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; 1 INTRODUCTION 



Millisecond pulsars (MSPs) can be used for a variety of as- 
trophysical applications and tests of ph ysics. They are sen - 
sitive probes of the interstellar medium (Backer et al. 1993) 



J1045— 4509 and J1643— 1224 were discovered during the 
Parkes survey for milliseco nd pulsars and, with the excep- 
tion of PSR J0437-4715 ( [Bell et al. 1995| ), have had only 
initial timing parameters reported, based on a data span of 



They provide mlormation about tnc evolution ot binary sys- 



tems and while dwarfs' (Bhattachafya & van den Heave] 



appr oximately one year (Bailes et al. 1991; Lorimer et al 



1991), and have unsurpassed sensitivity when used to search 
for planetary mass bodies (Wolszczan 1994). As excellent 



clocks, they challenge t he stability o f terrestrial time stan- 
dards in the best cases ( Taylor 1991 ) , place limits upon the 
energy density of cosmological gravitational waves ( Etomani 
fe Taylor 198s| ; [Kaspi, Taylor fc Ryba 1994), and are use d 



in testing theories of relativistic gravity (Taylor et al. 1992). 

Millisecond pulsars are generally weak radio sources, 
and their short periods make them difficult to discover. Over 
the last five years several surveys have increased from four to 
35, t he number of millisecond pulsars known in the Galactic 
disk (Camilo 1996), more than half of these new discoveries 
resulting from the Parkes southern sky survey (Manchester 
et al. 1996; A. G. Lyne et al. in prep.) . 

In this paper we report on the timing results ob- 
tained from observations of four millisecond pulsars with 
data spans of over 1000 d. PSRs J0437-4715, J0613-0200, 



1995 ). The observations and timing analysis are described 
in Section ^| while the timing residuals, pulsar parameters, 
and pulse profiles are presented in Section |^. Observations 
of the time variation of the projected semi-major axis of the 
orbit of PSR J0437-4715 is reported in Section |. In Section 
we summarise all measurements of orbital eccentricities for 
binary pulsars while, in Section ^ we use Lomb-Scargle spec- 
tra to assess the timing residuals for evidence of planetary 
mass companions. 



2 OBSERVATIONS AND TIMING ANALYSIS 

From the time of their discovery until 1996 March, at in- 
tervals averaging about three weeks, pulse arrival times for 
the four MSPs were obtained using the ATNF Parkes ra- 
dio telescope with cryogenic receivers at centre frequen- 
cies of 436, 660, 1520, 1940 and 2320 MHz. To date, only 
relatively small datasets have been obtained at 660, 1940 



© 0000 RAS 



2 J. F. Bell 



and 2320 MHz. Orthogonal linear polarizations were ob- 
served using a 2 x 256 x 0.125-MHz filter bank at 436 and 
660 MHz while both hands of circular polarization were ob- 
served with a 2 x 64 x 5.0 MHz filter bank at 1520, 1940 and 
2320 MHz. After detection, the signals from the two polar- 
izations were added, one-bit digitised at appropriate sam- 
pling intervals (0.3 ms at 436 and 660 MHz, 0.1ms at 1520, 
1940 and 2320 MHz, 0.08 ms for PSR J0437-4715) and writ- 
ten to magnetic tape. Off-line, the data were folded at the 
topocentric pulsar period to produce mean pulse profiles for 
each frequency channel, typically with integration times of 
90-180 s. These profiles were then transformed to the Fourier 
domain, phase-shifted appropriately to compensate for dis- 
persive delays across each observing band, transformed back 
to the time domain and summed to form a final mean profile 
at each frequency. 

Pulse arrival times for two of the pulsars, J0613— 0200 
and J1643— 1224, have also been obtained using the 76-m 
Lovell telescope at Jodrell Bank with cryogenic receivers at 
408, 606, and 1404 MHz. Both hands of circular polariza- 
tion were observed using a 2 x 64 x 0.125-MHz filter bank 
at 408 and 606 MHz and a 2 x 32 x 1.0-MHz filter bank at 
1404 MHz. After detection, the signals from the two polar- 
izations were added, filtered, digitised at appropriate sam- 
pling intervals, dedispersed in hardware before being folded 
on-line, and written to disk. 

For both data sets, a standard pulse template was fit- 
ted to the observed profiles at each frequency to determine 
the pulse times-of-arrival (TO As). The TOAs, weighted by 
their individual uncertainties determined in the fitti ng pro- 
cess, w ere analyse d with the TEMPO software package ( Taylor 
& Weisberg 1989), usin g the DE200 e phemeris of the Jet 
Propulsion Laboratory (standish 1982) and the Blandford 



& Teukolsky (1976) timing model for binary pulsars. Using 
the measured TOAs and an initial set of parameters describ- 
ing the pulsar system, TEMPO minimizes the sum of weighted 
squared timing residuals, the difference between observed 
and computed TOAs, yielding a set of improved pulsar pa- 
rameters and post-fit timing residuals. The dispersion mea- 
sures were determined by measuring the delays between the 
arrival times obtained at Parkes at two frequencies. We also 
fitted offsets between all other sets of data for each pulsar, 
to account for the different pulse shapes for each observing 
system and at different frequencies. 




49000 49200 49400 49600 49800 50000 
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Figure 1. Post-fit timing residuals for four MSPs plotted against 
MJD. 



3 TIMING RESIDUALS, PULSAR 

PARAMETERS, AND PULSE PROFILES 

The parameters obtained from fits to the data for each pul- 
sar are shown in Table [lj. The uncertainties quoted for these 
parameters are twice the formal values reported by TEMPO. 
All parameters shown in Table |l| were included in a global 
fit, except for P, which was fitted for separately. The dis- 
persion measure variation of PSR J 1045— 4509 was also ob- 
tained from a separate fit. For PSR J0437-4715, only 1520 
to 2320-MHz data were included in the global fit due to the 
much larger TOA uncertainties at lower frequencies arising 
from the greater dispersion smearing and scintillation. In the 
case of PSR J1045-4509, only the 1520-MHz data were in- 
cluded in the global fit due to the nonmonotonic dispersion 
measure variation with time. For the other two pulsars all 



available data were used. Post-fit timing residuals are shown 
for the four MSPs as a function of time in Fig. |lj All radio 
frequencies are plotted with the same symbol since there are 
no frequency-dependent trends in the residuals. 

For PSRs J0613-0200, J1045-4509 and J1643-1224 
the newly obtained celestial coordinates and binary param- 
eters, in particular the eccentricities, are considerably im- 
proved over previously published values and are discussed 
in more detail in Section ^. The uncertainties in the proper 
motion for PSR J0437— 4715 are improved by a factor of 
10 while useful limits are obtained on the proper motions 
of the other three pulsars. The velocities shown in Ta- 
ble |l| were obtained using distances calculated with the Tay- 
lor & Cordes (1993) free-electron-distribution model. For 
PSR J0437— 4715, the upper limit on the parallax of 6 mas 
constrains the distance to this pulsar to be greater than 
166 pc (la), 90 pc (2a), compared to the dispersion measure 
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distance of 140 pc with a 30 per cent uncertainty. Assuming 
that the intrinsic period derivative is zero, so that the mea- 
sured period derivative is entirely due to t he time-varying 



Doppler effect caused by the proper motion (Shklovskii 1970 



Camilo, Thorsett & Kulkarni 1994), gives an upper limit on 
the distance of 208 pc. The upper limits on P for the four 
pulsars provide no evidence for timing noise detectable with 
the current level of timing precision and data spans. 

In some cases, deviations of the post-fit timing residu- 
als from zero are somewhat larger than would be expected 
from the independently estimated TOA uncertainties. We 
believe that the TOA uncertainties are slightly underesti- 
mated, and we account for this in the final fitting procedure 
by multiplying the nominal TOA uncertainties by a small 
factor (typically approximately two) so as to produce a re- 
duced x 2 f° r t ne fit of unity. The original uncertainties are 
shown in the figures. 

In the case of PSR J0437— 4715 these deviations are 
very large (Fig. |l|), and we believe they have a systematic 
origin. From one observation of this pulsar to another the 
ratio of the main peak intensity to that of its outlying com- 
ponents appears to vary by up to 20 per cent, while the ratio 
of the component immediately following the main peak to 
the underlying emission varies by factors as large as two. 
We do not know the exact cause of these variations, but 
suspect that the observed changing pulse shapes may be 
caused by a relative gain variation in the two polarization 
channels and the high degre e of polarization of the emissio n 
observed across the profile ( Manchester fc Johnston 1995 ). 



To obtain the best set of profiles for timing, the profiles for 
different days were aligned and summed using the best avail- 
able ephemeris. Each profile was then normalised and sub- 
tracted from the sum. The resulting difference profiles were 
used to select a sample of profiles most similar in shape to 
the mean profile, and were subsequently used to obtain the 
timing solution presented. In making this selection, some 
time-span has been sacrificed but the timing accuracy has 
been significantly improved. Nevertheless, the timing preci- 
sion obtainable for PSR J0437— 4715 with these data is not 
yet sufficient to measure parameters such as parallax, or the 
apparent orbital perio d derivative caused by the large proper 
motion of the pulsar (Bell & Bailes 1996). 



Integrated pulse profiles for each of the pulsars are 
shown in Fig. H at four frequencies. In all cases the in- 
trinsic pulse widths are narrower than the instrumental 
resolution. At 660 MHz the resolution for all four pulsars 
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Figure 2. Integrated pulse profile for each pulsar at four frequen- 
cies. In each case the full period is shown and the instrumental 
resolution is indicated by a horizontal bar. 



however, somewhat higher than the 0.7/is predicted by the 
Taylor & Cordes (1993) model. 

We have observed dispersion measure variations for 
PSRs J1045-4509 and J1643-1224 and obtained an upper 
limit on the variations for PSR J0613— 0200 as shown in Ta- 
ble [lj. In the case of PSR J1045-4509, DM was fitted over 
the first half of the data span only, where the variation is 
monotonic. To illustrate the variations for PSRs J1643— 1224 
and J1045-4509, the best fit to the Parkes 1520-MHz data 
was made. For PSR J1045-4509 the DM was fitted around 
MJD 49500 and for PSR J1643-1224 the DM was fitted 
around MJD 49100. The residuals for the 436-MHz Parkes 
data where then calculated and plotted with the 1520-MHz 
data in Fig. ^. 

Backer et al. (1993) summarized observations of dis- 
persion measure variations for 13 pulsars and found that 
the rate of change of dispersion measure, ADM oc 
This is somewhat flatter than the linear depen- 



is determined by the sampling interval, as is tne resolu- 
tion of Xtm PSR J043 7 -4715 profiles at 436 and f 520 MHz. 

The resolution of the remaining profiles is dominated by 
dispersion smearing across individual channels. Hence, for 
these pulsars, the profiles observed at 1520 and 436 MHz 
are broader than the corresponding profiles at 2320 and 
660 MHz. The asymmetry and tail of the PSR J1643-1224 
profile at 436 MHz suggests some evidence for scattering 
along the line of sight to this pulsar. The sc attering time- 
scale t s (Taylor, Manchester & Lyne 1993) as measured 



VDM. 

dence observed for diffractive scattering (Cordes, Pidwerbet 



sky & Lovelace 1986). Dispersion measure variations sam- 



from the pulse profile in Fig. j2|, is approxim ately 240 /is at 
436 M Hz. After scalin g to 1 GHz by u~ AA ( fcordes, Weil 
berg & Boriakoff 1985), r s ,iGHz = 6/is and is comparable to 



pie electron-density fluctuations on length scales of 10 - 
10 15 cm, which are 10 4 times larger than the scale of the 
fluctuations which lead to diffractive scattering. This differ- 
ence in variation with dispersion measure suggests that these 
two effects may occur at different locations along the line of 
sight, or have different driving mechanisms. Our results are 
in general agreement with ADM oc %/ DM and therefore the 
conclusions drawn by Backer et al., although the scatter of 
the data about this relation is large. 



the scattering tim e-scales of other pulsars with simil ar dis- 
persion measures (Taylor, Manchester fc Lyne 1993). It is 
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Figure 3. Dispersion measure variations in two pulsars. Tim- 
ing residuals from Parkes data at 436 MHz (filled squares) and 
1520 MHz (circles), relative to a timing model fitted to the latter. 
Note the systematic variations in the difference between the two 
frequencies. 



4 SECULAR VARIATION OF ORBITAL 
INCLINATION 

It has been recently noted that the proper motion of a pulsar 
can lead to a secular variation of the projected semi-major 
axis of th e pulsar's orbit, a p sini, a s the orbital incl i nation i 



1996 



Sandhu 



changes ( Arzoumanian et al 
et al. 1)95). The contribution to the time derivative of the 



Kopeikin 1996 



projected semi-major axis is 

x — 1.54 x 10 -16 :e cot £(— /x a sin O + fj,g cos f2), 



(1) 



where x = a p sini/c is the projected orbital light travel time 
in seconds, Q is the longitude of the ascending node of the 
orbit, and /j, a and are the prop er motions in r ight ascen- 
sion and declination, in masyr -1 ( Kopeikin 1996[ ). Both the 
proper motions can be measured, while Q cannot, so that a 
measurement of x allows a constraint to be placed on i and 
hence the mass of the companion and on many aspects of 
the evolution of these systems. 

The types of system for which this method works well 
are those that have wide, small inclination (i.e. face-on) 
orbits, that are relatively nearby. It is therefore comple- 
mentary to the measurement of the "Shapiro delay" , which 
produces measurable effects for high-inclination (edge-on) 
systems in tight orbits. For millisecond pulsars with white- 
dwarf companions, the orbit needs to be very close to edge 
on for the Shapiro delay to be measurable. For intermedi- 
ate inclination angles, where some evidence of Shapiro de- 
lay is visible but where the covariance between the fitted 
companion mass and inclination angle is great (e.g. in the 
PSR J1713+0747 system, with i ~ 70°, Camilo, Foster, & 
Wolszczan 1994), measurement of x may be possible and 
should provide an additional constraint on the system. 



From the observations of PSR J0437— 4715, the mea- 
sured value of x — (7.4 ± 1.4) x 10~ 14 implies that x > 
Due to the non-Gaussian distribution of the tim- 



6x 10" 

ing r esiduals, we used bootstrap resampling (Press et al 



1992 ) to derive the uncertainties on all parameters in Table 



[]J including x. For most parameters, this yields similar un- 
certainties to those calculated in the standard way from a 
least-squares fit. However, for PSR J0437— 4715, the uncer- 
tainties on x and x are eight times larger using bootstrap 
resampling. We therefore take the more conservative boot- 
strap uncertainties as the 68 per cent confidence level. The 
maximum contribution from the proper motion occurs when 
the proper motion is aligned with the semi-major axis of the 
orbit projected onto the plane of the sky. This yields a limit 
of coti > 6.5 x 10 1 x/(xfi) = 0.83, giving an upper limit of 
i < 50°. Consequently, the companion mass is constrained 
to be greater than 0.19M©, for an assumed pulsar mass of 
1.4 Mq, considerably more than the a priori minimum pos- 
sible value of 0.14 M©. 

One may easily show that the measured x cannot 
be due to a genuine shrinking of the orbit. Since x — 
dp/csinz + etp/ccosi di/dt, one must ensure that a p sini <C 
a p cosidi/dt. This can be checked by obtaining an observa- 
tional limit on Pb and differentiating Kepler's third law with 
respect to time, which gives a p as a function of Pb. Inserting 
the observed limit for PSR J0437-4715 of P b < 2.2 x 10" 11 
shows that x < 3.1 x 10 -15 which is indeed much smaller 
than the measured value of x, and hence shows that x ~ 
a p /c cos idi/dt. The negligible value of a p is exactly as ex- 
pected since orbital evolution due to gravitational radiation 
or tidal effects leads to values of a p many orders of magni- 
tude smaller than observed for x. 



5 ORBITAL ECCENTRICITIES 

The timing observations presented here, vastly improve the 
accuracy of the eccentricity measurements, particularly for 
PSRs J0613-0200 and J1643-1224. In this section, we col- 
late all such measurements of, and upper limits on, eccen- 
tricities of binary pulsar orbits and compare them with the- 
oretical predictions. 

The convective fluctuation-dissipation theory of Phin- 
ney (1992) predicts a strong correlation between the orbital 
eccentricities and orbital periods of binary pulsars that have 
been recycled by stable mass-transfer from a Roche-lobe- 
filling, low-mass red giant. The measured eccentricities to- 
gether with upper limits are shown for all Galactic binary 
pulsars in Fig. ^j. All measurements (circles) and upper limits 
(squares) for objects with appropriate evolutionary histories 
are in excellent agreement with the model. As shown by the 
dashed lines, the trend of the upper limits for those objects 
with short orbital periods is merely a reflection of the timing 
precision attainable. 

Neutron stars in highly eccentric orbits have provided 
some remarkab le tests of general relativity and other theo- 
ries of gravity (Taylor et al. 1992). There are also interest- 



ing and useful tests than can be done using neutron stars 
in ver y circular orbits: tests of the strong equivalence prin- 
ciple ( Damour fc Schaier 1991 ), o f local Lorentz invariance 
( Damour fc Esposito-Farese 1992 ) and of conservation laws 
(Bell fc Damour 1996). These latter two tests depend lin- 
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perturbations in the timing residuals (Wolszczan 1994). Nat 




1 2 

log 10 (Orbital Period) 

Figure 4. Orbital eccentricities from Table 2 of Camilo (1996) 
and this paper. Solid line — the predicted median eccentricity for 
pulsars recycled by stable mass transfer from a Roche-lobe- filling, 
low-mass red giant. Dotted lines — 95 per cent of the measured 
eccentricities for recycled pulsars arc expec ted to lie within these 
(Phinney 1992; Phinney & Kulkarni 1994). Circles — low-mass 



binary pulsars. Squares — upper limits on measured eccentricities. 
Crosses — high-mass binary pulsars. Triangles - intermediate-mass 
binary pulsars. Dashed lines - approximate limits obtainable on 
eccentricities for timing residuals of 10 fis and 1 fis. 



early on the orbital eccentricity of a sample of binary mil- 
lisecond pulsars. 

The orbit with the lowest eccentricity in Fig. |L aniL 



in fact the known universe, is PSR J2317+1439 (Camilo 



Taylor 1996). This has been used to obtain a 



Nice 

tight li mit on the local Lorentz invariance of gravity (Bell 



very 



Camilo & Damour 1996). However, due to an unfavourable 



orientation of its position and proper motion with respect 
to the dipolar structure of the cosmic microwave back- 
ground, the limit is not as tight as it might otherwise be. 
The predicted eccentricities for some millisecond pulsars 
(e.g. PSR J0613— 0200) are somewhat lower than that of 
PSR J2317+1439, suggesting that these systems may ulti- 
mately provide more stringent tests of local Lorentz invari- 
ance, and may also allow the tests of the strong equivalence 
principle and conservations laws to be improved. However, 
as Fig. ^ shows, this will require timing precisions of about 
1 [is, which will be very difficult to achieve. Nevertheless, 
any improvement on timing precision for the pulsars with 
only upper limits on eccentricity will provide further tests 
of the fluctuation-dissipation theory. 



6 LIMITS ON PLANETARY MASS 
COMPANIONS 

The discove ry of planetary mass co mpanions to a millisec- 
ond pulsar (Wolszczan & Frail 1992) was unexpected. Their 
existence has been confirmed by the detection of three-body 



urally, such a discovery stimulates interest in finding other 
MSPs with planetary mass companions. In this Section, we 
place limits on the existence of any planetary mass compan- 
ions to the pulsar discussed in this paper using Lomb-Scargle 
spectra. While these pulsars have stellar-mass companions 
and may be less likely to have planetary companions, we 
apply the same technique to four single pulsars in another 
paper (M. Bailes et al. in prep.). 

From the residuals shown in Fig. [l] there appears to 
be no evidence for modulations that could be attributed to 
planetary mass companions. The mass function for a pulsar 
binary is given by 



/ = 



(ra c sm i ) 
(m p + m c ) 2 



4ir 2 (a p sini) 3 



(2) 



where i is the inclination angle of the orbit, a p is the semi- 
major axis of the pulsar's orbit, Pb is the orbital period, G is 
the gravitational constant and m p and m c are the pulsar and 
companions masses. Assuming circular orbits and m c <C m p , 
rearranging equation |^ suggests upper limits on the mass of 
undetected companions of 



m c sm % — 



An 2 m 2 



0.0143M©AT a 



1/3 



lday 



-2/3 



(3) 



(4) 



where AT a is the timing signal (peak-peak delay) due to 
the light travel time across the pulsar's orbit and is the 
mass of the Earth. If no modulation is seen in the residuals, a 
limit on AT a of an unseen companion can be obtained from 
the variance of the residuals, allowing a limit to be placed 
on m c sini. For long orbital periods the sensitivity drops as 
P 6 3 due to the strong covariance between the orbit al period 



of a planet and the period derivative of the pulsar ( Thorsett 



Phillips 1992) 



The above limits should be considered as qualita- 
tive. Quantitative limits on the presence of modulations 
can be obtained through a number of approaches: fast 
Fouri er transforms (FFTs), fol d and average, fitting of or- 
bits (frhorsctt fc Phillips 1992|), and Lomb- Scargle spectra 
flLomb 197$ |Scargic 1982| ; press et al. 1992[ ). Each of these 
methods has various advantages and disadvantages which 
are now briefly discussed. Fold and averaging, and fitting 
orbits to residuals are computationally expensive, provide 
poor detectability and sensitivity and do not offer a statis- 
tical assessment of the significance of any signals detected. 
FFTs overcome all of these difficulties, but suffer from alias- 
ing and poor sensitivity to long orbits resulting from the 
need to fill in gaps in the unevenly sampled data. 

The Lomb-Scargle periodogram overcomes these diffi- 
culties and is similar in speed to the FFT when one accounts 
for the additional data points required to provide the even 
sampling for the FFT. It also allows one to search for fre- 
quencies up to several times the mean Nyquist frequency. 
However, both the Lomb-Scargle periodogram and the FFT 
fail to take account of the covariance with other effects such 
as the Earth's orbit, w hich causes a loss of sensit ivity for pe- 
riods around one year ( Thorsett fc Phillips 1992 ). A distinct 
advantage of the spectral methods is that they do allow a 
formal assessment of the statistical significance of spectral 
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features, taking account of the number of frequencies that 
were searched. An ideal analysis therefore might consist of a 
Lomb-Scargle periodogram, combined with a small selection 
of fitted orbits to account for the above covariance problems. 

In this paper, examples of Lomb-Scargle spectra of the 
timing residuals for eight pulsars are given in Figs. ^ and ^| 
Four of the pulsars are those discussed in this paper and the 
remaining four are discussed in Camilo et al. (in prep.). The 
spectra as a function of frequency were converted to a func- 
tion of period, while the timing signal was obtained using 
AT a = 2 x 2(P a /N ) - 5 , where P a is the spectral power and 
No is the number of independent data points (Scargle 1982). 
Equation ^ was then used to convert to units of m c sini. To 
account for the covariance with the period derivative of the 
pulsar, the limits were multiplied by P^ for periods greater 
than one year. 

To test the sensitivity of the method, the eccentricity 
of PSR J1643— 1224 was adjusted by 0.24 per cent which is 
precisely three times the uncertainty in the measured value 
shown in Table ^ The resulting Lomb-Scargle spectrum in 
the lower panel of Fig. ^ shows a highly significant detec- 
tion. The probability that this signal is due to random noise 
is 1.5 x 1CP 5 . As another test of this method, the Jodrell 
Bank timing residuals for PSR B1257+12 were analysed in 
the same manner. After fitting for a slow-down model of 
the pulsar, the spectrum of the residuals (bottom panel of 
Fig. h) shows highly significant signals at the periods of the 
planets. The probability that these spikes are due to ran- 
dom noise is 2.8 x 10 -8 . It is important to note that the side 
lobes resulting from the uneven sampling of the data are 
very strong at all periods. For example, in the B 1257+ 12 
spectrum, the rms is totally dominated by the side lobes, 
while the true noise level for this pulsar is similar to that of 
the other pulsars shown in Figs. |E| and ^| This slightly in- 
creased rms may also be seen by comparing the spectrum for 
PSR J1643-1224 with the fake spectrum below it in Fig. |[ 

There are no spectral features that give any hint that a 
detectable periodic signal is present in the residuals of any 
of the eight MSPs considered. Probabilities of 0.1 and 0.999, 
that the signals are due to random noise are shown by the 
smooth curves in the figures, indicating that to a very high 
level of significance, the timing residuals are purely random. 
For the eight MSPs discussed here, companions with masses 
and orbits similar to Venus, Earth and the two largest plan- 
ets around PSR B1257+12 are ruled out. For those MSPs 
with more precise timing results, such as PSR J0437— 4715, 
the existence of companions with masses and orbits similar 
to Mercury and Mars are also ruled out. 




1 2 
log 10 (P b ) (days) 



Figure 5. The top four panels show Lomb-Scargle spectra of 
timing residuals for pulsars in this paper. In each panel the two 
smooth curves indicate the 0.1 and 0.999 probabilities that a sig- 
nal is due to random noise. Thus features above the line are 
likely to be significant signals and those below the line, noise. 
The lower panel shows the result of modifying the eccentricity 
of PSR J1643— 1224 to introduce a fake sinusoidal signal at half 
the orbital period. A significant detection is seen at precisely this 
period. 
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Table 1. Parameters of millisecond pulsars J0437-4715, J0613-0200, J1045-4509, and J1643-1224. Figures in parentheses represent 
la uncertainties in the last significant digits quoted. Some values of To and u> are given to greater precision than the available accuracy 
to help observers. 





-Ton ju^to/ — 4/10 


r OJA JUDIO — UZUU 


pop Tlflzl^ A^flQ 

roil j lu^to — ^touy 


pcR T~l fiA^ 1 99A 


a (J2000) 




04 h 37 m 15? 73501(2) 




06 h 13 m 43?97317(14) 




10 h 45 m 50? 1940(3) 


— 

16 h 43 m 38?1548(2) 


8 (J2000) 


—47° 15'08'.' 1527(2) 


— 02°00'47'.'072(6) 


— 45°09'54'.'204(2) 


— 1 2° 24' 58'.' 740 (14) 


H a (masyr 1 ) 


121.3(2) 


1.5(30) 


-7(3) 


7(3) 


US (masyr ) 


-70.4(3) 


-4(6) 


8(2) 


7(17) 


7r (mas) 


1(5) 








P (ms) 


5.75745182325672(13) 


3.0618440360787(6) 


7.4742241050635(10) 


4.6216414453446(9) 


P 


5.7296(8) x lO" 20 


9.57(4) x 10~ 21 


1.740(8) x 10~ 20 


1.840(6) x 10~ 20 


|P|(s _1 ) 


< 1.4 x 10~ 30 


< 5 x 10~ 3() 


< 14 x 10~ 30 


< 8 x 10~ 3 ° 


Epoch (MJD) 


49615.0 


49512.0 


49597.0 


49524.0 


DM (cm' 3 pc) 


2.649(4) 


38.7911(5) 


58.1833(18) 


62.4133(4) 


DM (cm -3 pcyr -1 ) 




0.0003(4) 


-0.0049(13) 


0.0010(5) 


Pb (d) 


5.741042355(3) 


1.198512557(3) 


4.083529188(17) 


147.017395(2) 


x ( s ) 




J- .uy ± ^iftu ^ u j 


q ni o$i(A\ 

O.VIO IZO^J 






7 A(~\ 4\ y 1 o — 14 


90 v 1 0~ 14 


QiAA\ v 1 O - 14 
y^iij x ru 


1 Pt v 1 O - 14 


e 


1.86(3) x 10- 5 


7(10) x 10" 6 


2.4(3) x 10" 8 


5.057(4) X 10~ 4 


U! 


2? 1811 ± 0.72 


43? 0440 ±36.5 


226? 078715 ± 10 


321?8578 ±0.04 


T (MJD) 


49615.852918±0.01 


49512.409269±0.2 


49598.199375±0.1 


49577.972054±0.015 


Pb 


< 2.2 x 10" 11 








fi (masyr -1 ) 


140.2(3) 


4(6) 


11(3) 


10(13) 


d (kpc) 


0.14 


2.2 


3.2 


>4.9 


v (km s~ 1 ) 


93(30) 


40(60) 


168(45) 




RMS residual (/is) 


1.1 


15 


10 


14 


Frequencies in fit (MHz) 


1520, 1940 & 2320 


400, 600 & 1400 


1400 


400, 600 & 1400 


Data span (MJD) 


49113-50117 


49030-50138 


49078-50117 


49046-50139 
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